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We present the design of a superconducting qubit that has 
circulating currents of opposite sign as its two states. The 
circuit consists of three nano-scale aluminum Josephson junc- 
tions connected in a superconducting loop and controlled by 
magnetic fields. The advantages of this qubit are that it can 
be made insensitive to background charges in the substrate, 
the flux in the two states can be detected with a SQUID, and 
the states can be manipulated with magnetic fields. Coupled 
systems of qubits are also discussed as well as sources of de- 
coherence. 



I. INTRODUCTION 

Quantum computers are devices that store information on 
quantum variables such as spins, photons, and atoms, and 
that process that information by making those variables in- 
teract in a way that preserves quantum coherence [Q-^] . Typ- 
ically, these variables consist of two-state quantum systems 
called quantum bits or 'qubits' Q. To perform a quantum 
computation, one must be able to prepare qubits in a de- 
sired initial state, coherently manipulate superpositions of 
a qubit's two states, couple qubits together, measure their 
state, and keep them relatively free from interactions that in- 
duce noise and decoherence ■ Qubits have been phys- 
ically implemented in a variety of systems, including cavity 
quantum electrodynamics [g| , ion traps , and nuclear spins 
Essentially any two-state quantum system that can 
be addressed, controlled, measured, coupled to its neighbors 
and decoupled from the environment, is potentially useful for 
quantum computation and quantum communications [ |l3| , [L4[ . 
Electrical systems which can be produced by modern lithogra- 
phy, such as nano-scaled quantum dots and tunnel junctions, 
are attractive candidates for constructing qubits: a wide va- 
riety of potential designs for qubits and their couplings are 
available, and the qubits are easily scaled to large arrays which 
can be integrated in electronic circuits [[],[l5| . For this reason 
mesoscopic superconducting circuits of ultra-small Josephson 
junctions have been proposed as qubits Jl6|-po[ and we detail 
one such circuit in this paper. 

Compared with the photonic, atomic, and nuclear qubits 
already constructed, solid state proposals based on lithog- 
raphy such as the one described here have two considerable 
disadvantages and one considerable advantage. The first dis- 
advantage is noise and decoherence the solid state en- 



vironment has a higher density of states and is typically more 
strongly coupled to the degrees of freedom that make up the 
qubit than is the environment for photons in cavities, ions in 
ion traps, and nuclear spins in a molecule or crystal. Extra 
care must be taken in solid state to decouple the qubit from 
all sources of noise and decoherence in its environment. The 
second disadvantage is manufacturing variability js| : each ion 
in an ion trap is identical by nature, while each lithographed 
Josephson junction in an integrated circuit will have slightly 
different properties. Solid state designs must either be insen- 
sitive to variations induced by the manufacturing process, or 
must include a calibration step in which the parameters of 
different subcircuits are measured and compensated for Jl5|. 

The advantage of solid state lithographed circuits is their 
flexibility: the layout of the circuit of Josephson junctions or 
quantum dots is determined by the designer, and its param- 
eters can be adjusted continuously over a wide range. As the 
results presented in this paper demonstrate, this flexibility al- 
lows the design of circuits in which the variables that register 
the qubits are only weakly coupled to their environment. In 
addition, the flexibility in circuit layout allows many possible 
options for coupling qubits together, and for calibrating and 
adjusting the qubits' parameters. That is, the advantage of 
flexibility in design can compensate for the disadvantages of 
decoherence and manufacturing variability. 

The flexibility in design afforded by lithography conveys a 
further advantage to constructing quantum computers. As 
noted above, a qubit has to accomplish at least five func- 
tions: it has to be addressed, controlled, measured, coupled 
to its neighbors and decoupled from the environment. One 
of the axioms of design, is that the number of parameters 
that characterize a system's design should be at least as great 
as the number of parameters that characterize the system's 
function |^l[ . The problem of having too few design param- 
eters available is particularly acute in the design of quantum 
computers and qubits : a quantum computer is a device in 
which essentially every physical degree of freedom is used to 
register information and to perform the computation. De- 
grees of freedom that are not used to compute are sources of 
noise and must be isolated from the computing degrees of free- 
dom. Designs for quantum computers are accordingly more 
constrained by fundamental physics than are designs for con- 
ventional computers: if one is storing information on a cesium 
atom, then the "design parameters" of the cesium atom — its 
energy levels, decoherence times, interaction strengths, etc. 
— are fixed by nature once and for all. In the lithographed 
Josephson junction circuits proposed here, by contrast, it is 
possible to make qubits that have a variety of different de- 
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sign parameters, each of which can be adjusted to optimize 
different functions. 



II. JOSEPHSON-JUNCTION QUBITS 

The superconducting Josephson tunnel junction is de- 
scribed by a critical current I and a capacitance C. (We 
will assume that the resistive channel of the junction is neg- 
ligibly small.) For superconducting circuits the geometrical 
loop inductance L s is also important if A = Lj/L 3 < 1, where 
Lj = & /2irI is the inductance associated with a Josephson 
junction in the loop. Here $„ = h/2e is the superconducting 
flux quantum. Josephson circuits can be divided into two gen- 
eral categories. Circuits of the first type have A > 1 so that 
the induced flux in the loop is not important. These circuits 
are typically made of aluminum, and the mesoscopic nature 
of their electronic transport has been studied in nanoscaled 
circuits. Circuits of the second type have A <C f and induced 
flux caused by circulating currents is important. These cir- 
cuits are typically made of niobium, and the macroscopic na- 
ture of the tunneling of flux has been studied in small-scaled 
circuits. 

The prospects of using superconducting circuits of the first 
type as qubits is encouraging because extensive experimental 
and theoretical work has already been done on mesoscopic su- 
perconducting circuits. (For a review of this work see Chap- 
ter 7 in J22j and in Ref. |^.) In circuits of the first type 
(A> I), two energy scales determine the quantum mechani- 
cal behavior: The Josephson coupling energy, Ej = I <fr /2ir, 
and the coulomb energy for single charges, E c = e 2 /2C. The 
energies can be determined by the phases of the Cooper pair 
wave function of the nodes (islands) and the number of excess 
Cooper pairs on each node. The phase and the number can 
be expressed as quantum mechanical conjugate variables 

In the "superconducting" limit Ej > E c , the phase is well 
defined and the charge fluctuates. In the "charging" limit, the 
charges on the nodes are well defined and the phase fluctuates 
strongly. When E.j and E c are within a few orders of mag- 
nitude of each other, the eigen states must be considered as 
quantum mechanical superpositions of either charge states or 
phase states. Such superposition states are important in de- 
signing qubits. Experimental studies have been performed by 
several groups with aluminum tunnel junctions with dimen- 
sions below fOOnm |2^^|. Superposition of charge states 
in circuits in the charging regime have been demonstrated 
[p5|-]27[ and are in quantitative agreement with theory pj| , p9[ . 
The Heisenberg uncertainty principle has been demonstrated 
when Ej w E c [ pOpo) . When Ej > E c topological excitations 
known as vortices exists and quantum mechanical interference 
of these quantities has been observed |n| . Unfortunately cir- 
cuits of the first type in the charging regime are sensitive to 
fluctuating off-set charges that are present in the substrate 
[ p2|j3^ ]. These random off-set charges make difficult the de- 
sign of a controllable array of quantum circuits and introduce 
a strong source of decoherence. 

In circuits of the second type (A <C 1), the quantum vari- 
ables can be related to the flux in the loops and their time 
derivatives. For a superconducting loop with a single Joseph- 
son junction, known as an rf SQUID, thermal activation of 
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FIG. 1. The three-junction qubit. Josephson junctions 1 
and 2 both have Josephson energies Ej and capacitance C 
and Josephson junction 3 has a Josephson energy and ca- 
pacitance a times larger. The nodes 1 and 2 represent the 
superconducting islands (nodes) which are coupled by gate 
capacitors C g = jC to gate voltages Va and Vb- The arrows 
define the direction of the currents. The flux is taken out of 
the page. 



macroscopic quantum states |34|] has been observed as well as 
macroscopic quantum tunneling between states and the dis- 
crete nature of the quantum states J35| . One of the advantages 
of these rf SQUID systems is that the two states have circu- 
lating currents of opposite sign and, hence, produce a readily 
measurable flux of opposite signs. To date no superposition 
of states have been demonstrated in the niobium circuits, al- 
though the improving quality of the niobium tunnel junctions 
may allow such a demonstration | p6|j37| . 

The goal of this paper is to design a qubit using circuits 
of the first type with aluminum, yet to have states (like in 
circuits of the second type) that are circulating currents of 
opposite sign. These circulating current states create a mag- 
netic flux of about 10 _3 $o and we refer to these as "persistent 
current (PC) states." These states obey all five functional 
requirements for a quantum bit: (1) The superconducting 
circuit is at sufficiently a low temperature that the the PC 
states can be made insensitive to background charges and ef- 
fectively decoupled from their electrostatic environment. The 
magnetic coupling to the PC states and the environment can 
also be made sufficiently weak. 



III. THE CIRCUIT 

The circuit of the qubit is shown in Fig . |l| Each junction 
is marked by an "x" and is modeled J22[|3q ] by a parallel com- 
bination of an ideal Josephson junction and a capacitor d. 
The parallel resistive channel is assumed negligible. The ideal 
Josephson junction has a current-phase relation, Ii = I sin y>i 
where ipi is the gauge-invariant phase of junction i. 

For the calculation of the energy the inductance of the loop 
is considered negligible A> 1 so that the total flux is the ex- 
ternal flux. In this case, fluxoid quantization around the loop 
containing the junctions, gives tpi — y>2 + fa — —2nf. Here 
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FIG. 2. U /E.j vs. / for a — 0.8 and for minimum energy 
phase configuration. The energy is periodic with period / = 1 
and is symmetric about / = 1/2. Near / = 1/2, there is a 
region [1/2 — f c , 1/2+ f c ] where there are two stable solutions. 
The inset plots f c as a function of a. 



f is the magnetic frustration and is the amount of external 
magnetic flux in the loop in units of the flux quantum $ . 

The Josephson energy due to each junction is Ej n (l — 
cos<p n )- The total Josephson energy U is then U — 
y\ Eji (1 — cos <pi). Combined with the flux quantization con- 
dition the Josephson energy is |39| 
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2 + a — cos (pi — cos if2 ~ ol cos(27r/ + (pi — if 2) ■ (1) 



The important feature of this Josephson energy is that it 
is a function of two phases [^| . For a range of magnetic frus- 
tration /, these two phases, (pi and 952, permit two stable 
configurations which correspond to dc currents flowing in op- 
posite directions. We illustrate this in Fig. ^, where we plot 
the energy of the minimum of the system as a function of / 
for a = 0.8. The energy is periodic with period / = 1 and is 
symmetric about / = 1/2. Near / = 1/2, there is a region 
[1/2 — f a , 1/2 + f c ] where there are two stable solutions. The 
inset plots f c as a function of a. These two solutions have cir- 
culating currents of opposite direction and are degenerate at 
/ = 1/2. The calculation of the energy for the stable solutions 
and f c is given in Appendix A. 

The main feature of the qubit that is proposed in this paper 
is to use these two states of opposite current as the two states 
of the qubit. By adding the charging energy (the capacitive 
energy) of the junctions and considering the circuit quantum 
mechanically, we can adjust the parameters of the circuit so 
that the two lowest states of the system near / = 1/2 will 
correspond to these two classical states of opposite circulating 
currents. Moreover, we will show that these two states can 
be made insensitive to the gate voltages and the random off- 
set charges. The quantum mechanics of the circuit will be 
considered in detail in the next section. 

The stable classical solutions correspond to energy minima 
in U(ipi, (£2). Let's consider the case of / = 1/2. For a < 1/2, 
U has only one minimum at (pi — if 2 — mod 2-7T. Above 
the critical value of a — 1/2, this minimum bifurcates into 
two degenerate minima at fi = —f2 = if* mod 2-7T where 
cos ip* = 1/2q. The minima form a two-dimensional pattern 
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FIG. 3. (a) A contour plot of the Josephson energy (poten- 
tial energy) U((pi,(p%) for / = 1/2 for a = 0.8. The nested 
nearly circular shapes mark the maxima in the potential, and 
the figure-eight shaped contours enclose two minima, (b) A 
plot of the potential vs. if m , the phase along the direction 
between these two minimum in the same unit cell, (c) A plot 
of the potential vs. (p n , the phase along direction from one 
minima to its next nearest neighbor. Note that the barrier is 
a saddle point. The upper curve in each figure is for a = 1.0 
and the lower for a = 0.8. 



with the two minima at (if* , —if*) and (—if* , if*) repeated in 
a two-dimensional square lattice. This pattern can be seen 
in Fig. ^| which is a contour plot of the Josephson energy as 
a function of the phase variables for a — 0.8. The nested 
nearly circular contours mark the maxima in the potential. 
The figure-eight shaped contour encloses two minima. Fig. |^b 
shows the potential along f m , the phase between the two 
minimum in a unit cell; that is, along the line if 2 = —(pi. The 
upper curve is for a — 1.0 and the lower for a — 0.8. Fig. |^c 
shows the potential vs ip n , which connects one minimum (say 
at (—if* , ip*)) to its next nearest neighbor (at (ip*, 2n — if*)). 
For q = 0.8 the energy barrier between the two minima is 
much lower than the energy barrier from the minimum in 
one unit cell to the neighboring unit cell. For a = 1.0 the 
energy barrier from unit cell to unit cell is nearly the same 
as the barrier within the unit cell. The ability to manipulate 
the potential landscape by changing a will be important in 
designing the qubit. 

We now consider the electric energy T stored in the five 
capacitors in the circuit. Each capacitor Cj has a voltage 
across it of V* so that 



QgAVA — QgBVs ■ 



(2) 



Here j = 1, 2, 3 and gA and gB. The last two terms subtract 
the work done by the voltage source to give the available 
electric (free) energy H]. The voltage across each Joseph- 
son junction is given by the Josephson voltage-phase relation 
Vn = (<l > o/27r)(/p n , where the over-dot indicates a partial time 
derivative. The ground in the circuit labels the zero of poten- 
tial and is a virtual ground. 
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The voltage across the gate capacitor gA is V 3 a =Va — Vi 
and similarly for V g B = Vb — Vi. The electric energy can then 
be written in terms of the time derivatives of the phases as 



T = 



The constant term 
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C g ■ Vg has been neglected and 
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C g = jC 
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(4) 



(5) 



The classical equations of motion can be found from the 
Lagrangian C — T — U . We take the electrical energy as the 
kinetic energy and the Josephson energy as the potential en- 
ergy [Q.The canonical momenta is Pi = dC/dip\. To attach a 
more physical meaning to the canonical momentum, we shift 
the Lagrangian by a Galilean-like transformation to 



C = T-U 



(I) 



<P -Cg-Vg. 



The canonical momentum is then 
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and is directly proportional to the charges at the islands at 
nodes 1 and 2 in Fig. hi as 



(8) 



(For any Josephson circuit it can be shown that there exist lin- 
ear combinations of the phases across the junctions such that 
these linear combination can be associated with each node 
and the corresponding conjugate variable is proportional to 
the charge at that node taU44] . If self and mutual inductances 
are need to be included in the circuit (as we argue does not 
need to be done in our case), then additional conjugate pairs 
would needed Q.) 

The classical Hamiltonian, H = ^\ Piipi — C, is 



M 



+ U{0) (9) 



where the effective mass M = ($ /27r) 2 C is anisotropic and 
the induced charge on the island is Q 1 = C 9 ■ V g . When 
driven by an additional external current source, the classical 
dynamics of this system have been studied in recent years 
both theoretically |^,^] and experimentally p?|]48| 

Note that the kinetic energy part of this Hamiltonian is 
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(10) 



which is just the electrostatic energy written is terms of the 
charges and induced charges on the islands. Often this is the 
method used in discussing the charging part of the Hamil- 
tonian. See for example Reference |43| and the references 
therein. A characteristic charge is e and characteristic capac- 
itance is C so that the characteristic electric energy is the 
so-called charging energy, E c = e 2 /2C. 



IV. QUANTUM CIRCUIT 

The transition to treating the circuit quantum mechan- 
ically is to consider the classically conjugate variables in 
the classical Hamiltonian as quantum mechanical operators 
p9|j50| ]. For example, the momenta can be written as Pi = 
—ihd/dpi and P2 = ~ihd/d<p2 and the wave function can 
then be considered as >= 9(tpi,(p2). 

In this representation the plane-wave solutions, such as 
ip — exp{— i(£i(fi + £2^2)} corresponds to a state that has 
l\ Cooper pairs on island (node) 1 and 1% Cooper pairs on 
island 2. These plane-wave states are the so-called charging 
states of the system [ |5l|j2l| . Since a single measurement of 
the number of Cooper pairs on each island must be an integer, 
then so should the I's here. (Note the expectation value of 
the number of Cooper pairs is not restricted to an integer.) 
Furthermore, an eigen function ^(^1,(^2) can be written as 
a weighted linear combination of these charge states. This 
means that 9(tpi, tpz) is periodic when each of the phases are 
changed by 2n, as in the physical pendula |52Tj . 

By considering 9(<pi, 952) = exp{i(fci<pi + k^^jxifi, P2) 
with [ki,k 2 \ = -(7<7/2e)[Vk, Vb], the Hamiltonian for 
Xifit fz) is almost the same but the induced charges are now 
transformed out of the problem, and we refer to this new 
Hamiltonian as the transformed Hamiltonian Ht, where |53| 



H t = -P T -M~ 1 ■P + Ej{2 + a 
— acos(2-7r/ + (fii — 932)} ■ 



COS ipi — COS tf2 



(11) 



The resulting equation Htxifi, f>2) = E\(ip\, ip^) is the same 
as for an anisotropic, two-dimensional particle in the peri- 
odic potential U. The solutions are Bloch waves with the 
"crystal momentum" k-values corresponding to — k', which is 
proportional to the applied voltages. This choice of crystal 
momentum insures that ^{tpi, 952) is periodic in the phases. 

We will first present the numerical results of the energy 
levels and wave functions for the circuit. Then we will use 
the tight-binding-like approximation to understand the re- 
sults semi-quantitatively. 

The eigen values and eigen wave functions for the trans- 
formed Hamiltonian H t are solved numerically by expanding 
the wave functions in terms of states of constant charge or 
states of constant phase. The states of constant charge re- 
sult in the standard central equation for Bloch functions and 
are computationally efficient when E c > Ej. The states of 
constant phase are solved by putting the phases on a discrete 
lattice and the numerics are more efficient when Ej > E c . 
Since the Josephson energy dominates, we will show results 
computed using the constant phase states. (However, when 
we used the constant charge states, we obtained the same 
results.) 

The numerical calculations are done in a rotated coordi- 
nate system which diagonalizes the capacitance matrix C by 
choosing as coordinates the sum and difference of the phases, 
ip p — ((pi + if2)/2 and ip m = (tpi — V?2)/2. The resulting 
reduced Hamiltonian is 



IP 2 IP 2 

—2 cos ip p cos ip m — a cos(27r/ + 2<p m )} 



(12) 
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f=0.5, k p =0 



FIG. 4. The energy levels E vs. frustration and 
gate voltage for Ej/E c — 80, a — 0.8, and 
7 = 0.02. The gate voltage is related to the k values by 
[k p ,k m ] = (jC/2e)[V A + V B ,V A - V b ], (a) E/Ej vs f b near 
ft = 1/2 for [kp, k m ] = [0, 0], and (b) E/Ej vs k m for fc p = 0. 



where the momenta can be written as P p — —ihd/d<p p and 
P m = -iTid/dtpm- The mass terms are M p = (§ /2-k) 2 2C{\ + 
7) and M m = (<l> /27r) 2 2C(l + 2a + 7). In this coordinate 
system the full wave function ^/(ip p ,ip m ) = exp{i(k' p ip p + 

Kn'Pm)}x(<Pp> < Prn.) witn [ fc p, fc m] = -(jC/2e)[V A + V B ,V A ~ 

Vb] and H t x( i Pp, l Pm) = Ex{f P ,^p m )- Also the two minima 
of the potential U(ip p , <p m ) within a unit cell form a periodic 
two-dimensional centered cubic lattice with lattice constants 
ai = 27ria; and a2 = tti x + 7ri y . 

Fig. ^ shows the energy levels as a function of / and as 
a function of the gate voltage which is given in terms of k. 
We have taken Ej/E c = 80, a = 0.8, and 7 = 0.02 in this 
example. The energy levels are symmetric about / = 1/2. In 
Fig. ^a, we see that the two lowest energy levels near / = 1/2, 
have opposite slopes, indicating that the circulating currents 
are of opposite sign. We also see that there is only a small 
range of 0.485 < / < 0.5 where the qubit can be operated be- 
tween these states of opposite circulating current. This range 
is consistent with the range [| ± f c ] from the classical sta- 
bility as shown in Fig. |j| At / = 0.49 direct calculation of 
the average circulating current, < ^\I sin ipi \ty > gives that 
the circulating current for the lower level is I\/I = —0.70 
and for the upper level is I2/I0 ~ +0.70. (A calculation 
of the circulating current from the thermodynamic relation 
— <&o 1 9E n /df gives the same result.) For a loop of diameter 
of d — 10 /im, the loop inductance is of the order /j, d w 10 pH 
[§4). For I a » 400 nA ( Ej = 200 GHz ), the flux due to the 
circulating current is LI\ ~ 10~ 3 <& o , which is detectable by 
an external SQUID. Nevertheless, the induced flux is small 
enough, that we are justified in neglecting its effect when cal- 
culating the energy levels. 

The difference in energy between the lower and upper level 
at the operating point of / = 0.485 is about 0.1 Ej ~ 20 GHz. 
Moreover, Fig. §b shows that the energies of these levels is 
very insensitive to the gate voltages, or equivalently, to the 
random off-set charges. The numerical results show that the 
bands are flat to better than one part in a thousand, especially 
at / = 0.48. To understand the underlying physics, a tight- 
binding model is developed. 




FIG. 5. The eigen wave functions for the lower (^1) and 
upper (^2) energy levels at / = 1/2 as a function of the 
phases. 



A. Tight-Binding Model 

Consider the case near the degeneracy point / = 1/2. The 
minima in energy occurs when ip p = and ip m = ±y? m where 
cos = 1/2q. Near the minimum at [tfim, <fip] = [<Pm, 0], the 
potential looks like a double potential well repeated at lattice 
points ai = 27ri^ and a2 = Tri^ + it\ y . Fig. |^ shows the two 
eigen functions in a unit cell. The wave function for the lower 
level (^1) is symmetric and the wave function for the upper 
level (^2) is antisymmetric. Both of the wave functions are 
localized near the two minima in U in the unit cell. 

To find an approximate tight-binding solution, let 
u((p m ,<p p ) be the wave function for the ground state on one 
side of the double potential wells, and v(ip m , <p P ) be the wave 
function on the other side. The tight-binding solution for Ht 
in Eqn. [12l is $ = c u u + c v v and satisfies 



E 



(13) 



Because the double well is symmetric at / = 1/2, each wave 
function has the same energy eo and so H uu = H vv = eo • Let 
ti be the tunneling matrix element between these two minima 
in the same unit cell and ti between nearest neighbor minima 
in the adjacent unit cells. Then H uv — H* u — —ti~t2e lk ' R2 — 
t 2 e ik (ai " a2) . The eigen energy levels are E = e ^\H uv \. The 
effect of ti is to split the degeneracy of the two states so that 
at k = 0, the energy is eo T (2ta +ti) for the symmetric and 
antisymmetric states respectively. The effect of t% is to give 
dispersion in k, that is, in gate voltage and off-set charges, 
to the energy levels. Because we want to minimize the gate- 
voltage (and off-set charge) dependence, we seek to minimize 
the tunneling t2 from one unit cell to another. Likewise, we 
want the two localized states in the two wells to interact, 
so that we want t\ to be non-zero. This is why the potential 
landscape in Fig. ^| was chosen to have a w 0.8: The potential 
has a much lower barrier between states in the double well, 
but a large barrier between states from one double well to the 
next. 

An estimate of ti can be obtained from calculating 
the action Si between the two minima and using ti w 
{hi0i/2Tv)e- Si/n where LUi is the attempt frequency of escape 
in the potential well. The action from point (p a to (pi is 
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S = 



(2M nn (E - U)) 1/2 \dip n 



(14) 



Here n is a unit vector along the path ol integration, dip n 
the differential path length, and M nn — n T • M • n is the 
component of the mass tensor along the path direction. In 
both cases we will approximate the energy difference E — U 
as the difference in the potential energy AU from the minima 
along the path. 

First, consider the calculation of ti, the tunneling matrix 
element within the unit cell. The path of integration is taken 
from (— <pj^,0) to (ip m ,0) along the direction n = i x , so that 
M, ln = M m for this path. The potential energy at the minima 
is Umin = 2 — l/2a. The difference in the potential energy 
from the minima at (— ip m , 0) along this path is can we written 
as AUi = Ej{2a (cos f> m — 1/2q) 2 }. The action along this 
path is then 



Si 



(4M m aEj) 1/2 (costpm - 



(15) 



which yields 

Si = h [4q(1 + 2a + ~/)Ej/E c ] 1/2 (sin tp* m . (16) 

Now consider t 2 , the tunneling from unit cell to unit cell. 
For example, take the integration to be from (<p m ,0) to one 
of its nearest neighbor minima at (n — ip m ,n). We will take 
the path of integration to be a straight line joining these two 
points in the ifm-^fp plane. This path is not the optimal 
trajectory, but the difference of this straight line path from 
the optimal trajectory is quadratic in the small deviations 
of these two paths. The straight line path is described by 
<Pm = tp m + A<^)p where A = (it — 2(p m )/ir; it has a direction 
of n = Aia, + i y and a path length of ds = \/l + A 2 dtp v . 
The mass on this direction is M 2 = {M p + A 2 M m )/(l + A 2 ). 
The difference of the potential energy along this path form 
the minima energy is AU2/EJ — —2 cos ip p cos + Aip p ) + 
2a cos 2 ((p m + \<p P ) + l/2a. The action is then 



S 2 = [2M 2 Ej(l + \ 2 )] 1/2 

Jo 



( AU 2 
\ Ej 



1/2 



dip p . 



(17) 



The integrand is not analytically integrable, but it is zero at 
the end points of the integration and is well approximated by 
y/AU 2 /Ej w (1/V2a) cos(ip p - ir/2). With this approxima- 
tion, S 2 = (4M 2 Ej(l + A 2 )/a) 1/2 , which is 



S 2 



(18) 



To compare the tunneling rates we would first need the 
attempt frequencies in the two directions. However, we can 
consider the attempt frequencies to be of the same order of 
magnitude and thus t 2 /ti ~ e ~ (S2 " Sl)/R . For a = 0.8, we 
find that Si/{h^fEj/E c ) w 0.6 and S 2 /(h ^Ej/E c ) w 1.4. 
For Ej/E c ~ 100, then t 2 /ti ~ 10 -4 < 1. We are therefore 
able to ignore t 2 , the tunneling from the unit cell to unit cell. 
This means that there is little dispersion in the energy levels 
with k and consequently, with the voltage or off-set charges. 
In fact, using the action one can show that for a smaller than 



about 0.85, ti > t 2 for Ej/E c ^ 80. Throughout the rest of 
the paper we will choose parameters so that the effects of t 2 
can be ignored. 

We now obtain an approximate solution for the energy lev- 
els and tunneling matrix elements by modeling each side of the 
double potential. Near the minimum at [ip m ,ipp\ = [^^,0], 
the potential looks like an anisotropic two-dimensional har- 
monic oscillator. The Hamiltonian in the vicinity of the min- 
imum is approximately, (with Q\ = 0) 



1 Pi 



1 1 P 2 

H ^2M- + 2 M ^ + 2M m 



where 



and 



~e7 



a(l + j)(Ej/E c ) 



4(4a 2 - 1) 
*(l + 2a + >y)(Ej/E c ) 



+ U 



(19) 



(20) 



(21) 



and Uo — 2 — l/2a. The ground state <j> of the single har- 
monic well has energy eo = h(ujp + io m )/2 + U - Let's now 
use this approximation to understand the energy levels, first 
at / = 1/2 and then near this point. 

At / = 1/2 we expect the four lowest energy levels of the 
two-dimensional harmonic oscillator to be with u m < uj v , 
Ei — eo — ti, E 2 = eo + ti, E3 = eo — ti + hiv m , E4 = 
eo + ti + fnum ■ Table [j] compares the results and we also list 
the small anharmonic corrections to the simple harmonic en- 
ergy levels. We have chosen to compare (Ei + E 2 )/2 and 
(£3 + Ei)/2 so that the tunneling term is absent and a direct 
comparison with the simple harmonic oscillators can be made. 
The agreement between this tight-binding approximation and 







hujp 


Eo 


(£1 + E 2 )/2 


(E 3 + E 4 )/2 


Harmonic 


0.193 


0.247 


1.60 


1.79 


1.84 


Anharmonic 


0.183 


0.238 


1.59 


1.77 


1.83 


Numerical 


0.154 


0.226 


1.58 


1.74 


1.81 



TABLE I. A comparison of the energy levels with the 
approximate harmonic oscillator levels (with harmonic and 
anharmonic terms) with the numerical calculations. Here, 
/ = 1/2, a = 0.8, 7 = 0.02, and and Ej/E c = 80. Also 
Uo ~ 1-38 and [/bar = 0.225 for the harmonic and anharmonic 
estimations. All the energies are in units of Ej. 

the numerical calculations is good. We have also included the 
barrier height from one minimum to the other one in the same 
unit cell. 

If we estimate the attempt frequency for t\ as uj m , then 
we find that for the parameters in Table |. that the action 
calculation gives ti = lQ~ i E. J . From the full wave functions, 
we estimate t\ = (E 2 — E\)j2 w W~ 3 Ej. This discrepancy 
can be made smaller by noting that in the calculation of the 
action, we could more accurately integrate from the classical 
turning points in the potential rather than from the minima 
[p5|. However, for our purposes, the action expression will be 
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sufficient for quaiitative discussions, and we wiii use the fuii 
numerical calculations when estimating actual numbers. 

So far we have estimated the energy levels and tunneling 
matrix elements when / = 1/2. As / is decreased from / = 
1/2 the potential U changes such that one well becomes higher 
than the other, and the barrier height also changes. For the 
qubit we are mainly interested in the lowest two energy states 
of the system, so we now estimate the terms in tight-binding 
expression of Eqn. By defining the zero of energy as the 
average of the two lowest energy states at / = 1/2, we find 
that the Hamiltonian for these two states is 



H 



F -t 
-t -F 



(22) 



Here F is the energy change of each of the wells measured with 
respect to the energy of the wells at the degeneracy point; that 
is, F = (dU/df)Sf where U is the potential energy energy. 
Note that since we will be operating the qubit just below the 
degeneracy point / = 1/2, then F < 0. Also, t = ti + At, 
where t\ is the intracell tunneling matrix element calculated 
at the degeneracy point and At is the change. The eigen val- 
ues are Ai,2 = tV F 2 + t 2 where we have explicitly assumed 
that F is negative and t is positive. 

The eigen vectors are given as the columns in the rotation 
matrix 



D{9) 



cos (9/2 -sin 0/2 
sin 0/2 cos 0/2 



(23) 



where = — arctant/f\ For example, at the degeneracy 
point, F — 0, so that E = =pt and the eigen vectors are 
(1A/2, 1/^2) T and (-I/a/2, 1/^2) t . These are just symmet- 
ric and antisymmetric combinations of the single well wave 
functions, as expected. For / slightly below 1/2, we have 
\F\ > t, so « 0, and the energies are E = tVF 2 + t 2 « ±F. 
The eigen vectors are approximately (1,0) T and (0, 1) T , so 
that the eigen states are nearly localized in each well. 

It is more convenient to discuss the Hamiltonian and eigen 
states in the rotated coordinate system such that Hn = 
D T (8)HD(8). In the rotated coordinate system, the Hamil- 
tonian is diagonal with 



H D = - v / F 2 +t 2 



(24) 



and the eigen energies are E — ±V F 2 + t 2 and the eigen 
states are then simply spin down |0 >= (1,0) T and spin up 
1 >= (0, 1) T vectors. In other words, no matter what the 
operating field is, we can always go to a diagonal represen- 
tation; but the rotation matrix must be used to relate the 
simple spin up and down vectors to the linear combinations 
of the wave functions in the well. 



V. MANIPULATION OF THE QUBIT 

As noted above, the flexibility of the design of Josephson 
junction circuits affords a variety of methods for manipulating 
and controlling the state of qubits. In this section we show 
how the basic qubit circuit can be modified to allow precise 
control of its quantum states. To manipulate the states of 
the qubit, we need control over the properties of the qubit. 



For example, control over /, the magnetic field, allows one 
to change the operating point and F, the value of the energy 
difference between the two states. Control over the potential 
barrier height allows changing of the tunneling through t. For 
example, if the operating point of F and t are changed by 
AF and At, then the Hamiltonian in the rotated coordinate 
system is 



-y/FS + tl a z + AH D 



(25) 



where with O = — arctani /F , 
AHd = A_F(cos g z — sin a x ) — Arj(sin0 a z + cos a x ) . 

(26) 

The control over F can be done by changing /. The control 
over t can be done by changing the barrier heights. To control 
the barrier heights by external parameters, we replace the 
third junction by a SQUID which acts like a variable strength 
junction. The modified circuit of the qubit is shown in Fig. ^. 



v„ A c 




v„„ c„ 



J >< Cj E ]2 



FIG. 6. The four-junction qubit. Two junctions form a 
SQUID loop and have Josephson Energies and capacitance /3 
times larger than the other junctions 1 and 2 which both have 
Josephson energies Ej and capacitance C. The nodes A and 
B represent the superconducting islands which are coupled by 
gate capacitors C g = 7C to gate voltages Va and Vb- The 
arrows define the direction of the currents. The flux is out of 
the page. 

We now analyze this circuit since it will be used in all sub- 
sequent discussion of the qubit. Flux quantization around 
each of the two loops, gives (pi — (p2 + <fi3 = — 2tt/i and 
(^4 — ^3 = — 2-71-/2. The Josephson energy due to each junction 
is Ej n (1 cos tp n ). The total Josephson energy U is then 



— - =2 + 2/3 — 2 cos ip p cos <p„ 
Ej 



2(3 cos(nf a ) cos(2nf b + 2<p m ) 
(27) 



where tp p = (</?i + <f2)/2 and (p m = (tpi — ¥ 3 2)/2, and also 
fa = /2 and fb — fi + jil2. Hence we see that 2/3 cos(n f a ) 
plays the role of a in the three junction qubit, but now this 
term can be changed by changing f a = f%, the flux in the top 
SQUID loop. Likewise, /j, = fx + f%/2 plays the role of / in 
the three-junction qubit. The reduced Hamiltonian is then 
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2W P + 2it; +Ej{2+w 



-2 cos ip p COS (fim — 2(3 cos(7r f a ) cos(2n f b + 2ip m )} (28) 



where M p = ($ /27r) 2 2C(l + 7) and M m = ($ /27r) 2 2C*(l + 
4/3 + 7). 

To manipulate the parameters in the Hamiltonian let the 
magnetic fields change very slightly away from the some 
degeneracy point of /* and /I to a new operating point 
fi = fi + £ i an d f 2 = fi + £ 2 ■ Then J 1 changes from zero to 
F = nei + T2£2 and i changes to t = t\ + siei + S2£2, where 
n and s; are constants and t\ is the tunneling matrix element 
at the degeneracy point as found in the previous section. We 
take the operating point to be effectively in the regime where 
/ < 1/2 in Fig. [|, so that ei, 2 < 0. Hence, F Q < 0. Also, 
t is assumed to remain positive. In the new rotated frame 
with d = — arctant /F , the Hamiltonian given by Eqn. ^ 
is H D = -VF2 + tl<T z . 

Away from this new operating point, let fi ~ f° + 5i and 
/2 = /a +62- In the operation of the qubit, 5i <C £i| and Si 
usually will usually have a sinusoidal time dependence. Then 
F — F + riSi + r 2 5 2 and t = t + si5i + S2S2, so that AF = 
ri5i + r2fe and At — s\Si + S2^2- Then Hamiltonian in the 
rotated system with o = — arct an t /F is 



-J- = 2n y/1 - l/(4/3 2 ) . 



(32) 



For our example with (3 — 0.8, we have ri/Ej = 4.90. Es- 
timates obtained from the numerical calculations done by 
changing /1 and /2, give ri/Ej = 4.8 and r2/Ej = 2.4 in 
good agreement with Eqn. Bq in Appendix B. 

Likewise, from Appendix B we have that si = and S2 = 
Tjt\J Ej/Ec where r\ is of the order of unity. For the operating 
point we find r\ ~ 3.5. Therefore, changes in H due to changes 
in t\ go like a x . These tight-binding estimates for (3 = 0.8 give 
si = and S2/EJ — 0.03. Full numerical calculations for our 
example with si = and S2/EJ = 0.20. The agreement with 
the tight-binding results are good, although the tight-binding 
underestimates S2 with for these parameters. 

In summary, from the degeneracy point of /* = f% = 1/3, 
let the operating point be f° = /* + ei and /| = /j* + £2, 
so that F = r 1(6162/2) and t = ti + S2 62- Now consider 
the changes in field about the operating point such that fi = 
/* + Si and /2 = /| + 82- In the rotated frame where O = 
— arctan t /F , the Hamiltonian is 



H D = -VFj+tl a z + AH L 



(33) 



where 



H D = -VFf+¥oa z + AH D 



(29) 



where 



AH D = {riSx + r 2 <52)(cose o (T z — sin# a x ) 

— (si<5i + S2<52)(sin0 o u z + cos O cr x ) . (30) 

Hence we see that changes in the magnetic field from the 
operating point of f° and /| cause both a z and <j x types of 
interactions. 

To find the magnitude of these changes, we calculate the 
coefficients of change (n, T2, si and S2) most simply at the 
degeneracy point where ei = 0; that is, at the degeneracy 
point f° — f*. We choose the degeneracy point for the four- 
junction qubit at /* = 1/3 and /| = 1/3. This results in 
classically doubly degenerate levels. In fact, any choice which 
satisfies 2/* + /jf = 1 when the classical energy U has two 
minima, will also result in doubly degenerate levels. For ex- 
ample /1 = 1/2 and /o* = is also a possible and convenient 
choice. However, we prefer /* = /I = 1/3 for the following 
reason. The change in potential energy with f a gives 



dU 

dfa 



~2ir/3 sin irf a cos 2^?^ 



^Jf = -2tT 2 /3 COS 71^008 2^^ 



(31) 



The first order terms vanishes if /I = 0, resulting in the 
potential barrier always decreasing with changes in / 2 . On the 
other hand, if /I = 1/3, then the barrier height can be made 
to increase and decrease with changes in /2, thus allowing 
more control of the qubit. 

Now the coefficients of change (ri, V2, si and S2) can be 
estimated both from the numerical calculations and from the 
tight-binding model as shown in Appendix B. We find that at 
the degeneracy point of /1 = / 2 = 1/3, 



AHd — ri (81 + — ) (cos do <Jz — sin O cr x 
— S2<52 (sin 6> a z + cos0 o a x ) ■ 



(34) 



and n/Ej = - l/(4/3 2 ) and s 2 = T)t y/Ej/E c . 

A typical design for a qubit will have Ej/E c = 80, (3 = 0.8, 
7 = 0.02. We find from numerical calculations that t « 
0.005-Ej and 77 « 3.5, which agree well with our tight-binding 
estimates. We operate at f x =f 2 = 0.33 so that ei = £2 = 
— 1/300. (This is equivalent to operating the three-junction 
qubit at / = /1 +/ 2 /2 = 0.495 in Fig. ^.) Writing the energies 
as Ei — hvi , we have taken typical values of Ej — 200 GHz 
and E c — 2.5 GHz, and we find that t — 1 GHz and F = 
5 GHz (which gives a splitting between the two states of about 
10 GHz). The Hamiltonian is found to be 



— - = -0.025cr, + (4.0*i + 2.15 2 ) cr z 
Ej 



(0.4651+0.4152)^ 
(35) 



The numerical values used are from numerical calculations. 
These values agree well with the estimates used in Eqns. ^ 
and for the level splitting and the terms proportional to 
n; the terms proportional to S2 match to about 50%, due to 
the more sensitive nature of estimating the tunneling terms. 

The terms containing a x can be used to produce Rabi os- 
cillations between the two states by modulating 5i and 52 
with microwave pulses of the frequency of the level splitting 
of 2F = 10 GHz. One could arrange the values of 5i and 52 
to make the time- varying a 2 -term vanish. Then the operation 
of the qubit would be isomorphic to the NMR qubit. How- 
ever, our simulations show that such an arrangement couples 
higher energy levels and invalidates the simple two-state ap- 
proximation. This is due to the large matrix element between 
the ground state and the second excited state given by the 
change in potential due to varying 52. (It is interesting to 
note that similar coupling to higher levels occurs in qubits 
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FIG. 7. Coupling of qubits ^4 and B through the mutual 
inductance between (a) the lower regions of both, and (b) the 
lower region of A and the upper region of B. 



based on the RF SQUID and on simple charge states.) We 
propose to manipulate the qubit by varying Si which causes a 
Rabi oscillation through the a x term as well as a strong mod- 
ulation of the Larmor precession through the time varying <r z 
term. Because the Rabi frequency is much smaller than the 
Larmor frequency, the precession causes no problem for ma- 
nipulating the qubit. For S\ = 0.001 and 82 = 0, the Rabi 
frequency is about 90 MHz. We note that this mode of opera- 
tion is also possible with the three-junction qubit. Of course, 
it will not be possible to completely eliminate the deleterious 
effects of the 82 coupling, but the effect of this coupling can 
be greatly reduced if 82 is restricted below 0.0001. 

The varying magnetic fields 8\ and 82 can be applied lo- 
cally to the qubit by using a control line to inductively couple 
to the qubit. Moreover, if the the control line is driven by 
an Josephson oscillator, then the coupling circuit could be 
fabricated on the same chip. 



VI. INTERACTION BETWEEN QUBITS 

A variety of methods are available for coupling qubits to- 
gether. As noted in p^|Jl4| , essentially any interaction be- 
tween qubits, combined with the ability to manipulate qubits 
individually, suffices to construct a universal quantum logic 
gate. Here we present two methods for coupling qubits in- 
ductively as shown in Fig. The inductive coupling could 
either be permanent, or could be turned on and off at will by 
inserting Josephson junctions in the coupling loops. 

Fig. Ma shows one way of coupling two identical qubits. 
The lower portions of each qubit (the loops that contain the 
circulating currents) are inductively coupled. To a first ap- 
proximation we model the coupling as changing the flux in 



each of the two lower rings only through the mutual induc- 
tive coupling. (We ignore the the self-inductance which can 
easily be included.) The effective frustration in the lower 
loop of A, fi, is changed over the applied frustration /j 4 to 



of B is Ii . Similarly, 
Hamiltonian is 



Here the current in the lower loop 
fx = fx + Mlf/& . The coupled 



H A1 



h a (A a 



\fx B 



+ MI1I1 



(36) 



which is the sum of the Hamiltonians for each system plus a 
term due to the mutual inductive coupling. 

The inductively coupled contribution to the frustration is 
estimated to be of the order of 10 -3( 3> o which is much smaller 
than the applied frustration. Since each persistent current 
will inductively couple into the other qubit, this will produce 
changes in the Hamiltonian of the a z and a x type and these 
changes will be proportional to the sign of the circulating 
currents in the qubit. Hence, we expect the coupling to be 
described by an interaction Hamiltonian of the form, 



7-yint A I 

Bab = «i o" z ff« 



A B . . A B 

K20"z a x + +K3(T X a z 



(37) 



and a 

We have estimated magnitude of 



°z <? x types of coupling. 
Ki w O.OlEj. 

As Eqn. shows, the inductive coupling between the 
qubits can be made to be a substantial fraction of the qubit 
Larmor frequency. This is an attractive feature, as the cou- 
pling between two qubits sets the speed limit for how rapidly 
two qubit quantum logic operations can be performed in prin- 
ciple. In practice, it may be desirable to sacrifice speed of 
operation for enhanced accuracy: in this case, the inductive 
coupling could be designed to be smaller by decreasing the 
overlap of the inductive loops with the circuits. 

Coupling between qubits is similar to the coupling we en- 
vision between the qubit and the measurement circuits con- 
taining SQUID-like detectors. In its usual configuration, the 
SQUID is biased in the voltage state which produces a volt- 
age related to the flux through its detector loop. However, 
such a strong, continuous measurement on a qubit would de- 
stroy the superposition of states in the qubit and project out 
only one of the states. This problem can be circumvented 
by designing a SQUID such that it is current biased in the 
superconducting state and hence is not measuring the flux in 
its detector loop. When one needs to measure the qubit, the 
SQUID can be switched to its voltage state, for example, by 
applying a pulse of bias. The coupling from mutual induc- 
tance between the SQUID and the qubit will also have to be 
controlled. Other measurement schemes using SQUIDs which 
are weakly coupled to the macroscopically coherent system 
have been proposed [M. 



VII. COMPUTING WITH THE PC QUBIT 

All the ingredients for quantum computation are now avail- 
able. We have qubits that can be addressed, manipulated, 
coupled to each other, and read out. As will be indicated be- 
low, the particular qubits that we have chosen are well insu- 
lated from their environment as well. The flexibility of design 
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for collections of qubits now allows a wide variety of overall 
designs for quantum computers constructed from such qubits. 

Before discussing various superconducting quantum com- 
puter architectures, let us review some basic ideas about 
quantum logic and see how to implement quantum logic us- 
ing our superconducting qubits. A quantum logic gate is a 
unitary operation on one or more qubits. Quantum computa- 
tions are typically accomplished by building up quantum logic 
circuits out of many quantum logic gates. Just as in the case 
of classical computers, certain sets of quantum logic gates are 
universal in the sense that any quantum computation can be 
performed by wiring together members of the set. In fact, al- 
most any interaction between two or more qubits is universal 
[ p^|Jl^] ; but a convenient universal set of quantum logic gates 
widely used in the design of quantum algorithms consists of 
single qubit rotations and the quantum controlled-NOT gate, 
or CNOT §?[. 



So the basic quantum logic operations can be performed 
on our superconducting qubits in a straightforward fashion. 
Accordingly, it is possible in principle to wire groups of qubits 
together to construct a quantum computer. A variety of ar- 
chitectures for quantum computers exist, usually consisting 
of regular arrays of quantum systems that can be made to 
interact either with their neighbors or with a quantum "bus" 
such as a cavity photon field or a phonon field in an ion trap 
that communicates equally with all the systems in the array. 
Because of the flexibility inherent in laying out the integrated 
Josephson junction circuit, a wide variety of architectures are 
possible. A particularly simple architecture for a quantum 
computer can be based on the proposal of Lloyd jljg for ar- 
rays of quantum systems such as spins or quantum dots. 

C. Linear Chain of Qubits 



A. One-Qubit Rotation 

An arbitrary one qubit rotations can be written as e~ lcrt — 
cos t — i sin t a for some Pauli matrix a = aa x + ba y + ca z , 
where a 2 + b 2 + c 2 = 1. There are many ways of accomplish- 
ing a one qubit rotation: the ability to rotate the qubit by 
a precise amount around any two orthogonal axes suffices. 
Pursuing the analog with NMR, we choose a method that 
involves applying an oscillatory field applied at the qubit 's 
resonant frequency to rotate the qubit. 

The Hamiltonian for a single qubit (A) can be gotten from 
Eqn. |i[ Here we assume Ej = 200 GHz, 81 = 0.001 cos tut 
and 82 = 0, and the level splitting is ui — 10 GHz. Then, the 
Hamiltonian is 

H D (GHz) = 5 a z + 0.80 (cos tut) a z - 0.09 (cos cut) a x (38) 

The Rabi frequency is 90 MHz so that a 7r pulse would be 
about 20nsec. 

B. Two-Qubit Controlled NOT 

A controlled NOT is a two qubit quantum logic gate that 
flips the value of the second qubit iff the value of the first qubit 
is 1. That is, it takes: |00) -> 1 00) , |01) -> |01), 1 10) -> |11), 
and 1 11) — * 1 10). A controlled NOT can be combined with sin- 
gle qubit rotations to give arbitrary quantum logic operations. 
A controlled NOT can be straightforwardly implemented in 
the superconducting qubit system by exploiting the analogy 
with NMR. Suppose that two qubits A and B have been con- 
structed with an inductive coupling between their lower loops 
as in the first part of the previous section. Then the level 
splitting of qubit B depends on the state of qubit A, with 
values A_Eyt,o for A in the |0 > state and AEa,i for A in the 
|1 > state. When a resonant pulse corresponding of AEa,i/Ti 
is applied to qubit B, it will only change if qubit A is in its 
1 1 > state. Since the coupling between the qubits is consider- 
ably larger than the Rabi frequency, the amount of time that 
it takes to perform the controlled NOT operation is equal to 
the amount of time it takes to perform a tt rotation of a single 
qubit. 



Consider a linear array of qubits ABABABAB .... Let the 
bottom of each qubit be inductively coupled to the top of the 
neighbor to the left. Also let each type of qubit, A and B, 
have a slightly different Josephson energy. Each qubit also 
has the area of the top loop which is half that of the bottom 
loop. In the absence of the driving electromagnetic fluxes (the 
51), the Hamiltonian for the system can be generalized to be 
written as 

H = (WftCTfc + 2J k}k+1 alal +1 ) (39) 

k 

where huj k = ^ F 2 + t\ and J k ,k+i = «fc,&+i(»~i,k +fi,fc+i)/2. 
This problem then maps on the linear chain of nuclear spins 
which was shown by Lloyd || to be a universal quantum com- 
puter. The coupling needed to perform 7r/2 pulses is provided 
by the terms containing the Sj's. The nice feature of this lin- 
ear chain is that separate control lines for AC fields are not 
needed. The whole linear array can sit in a microwave cavity 
and be pulsed at the desired frequency. (The dc bias fields 
to ensure fa = f 2 = 1/3 will require at least two dc control 
lines). The frequencies needed are around 10 — 25 GHz with 
intervals of 1 GHz (and with resolution of about 0.1 GHz). We 
could make these numbers larger or smaller if needed. 

Details of computing with this are given in various refer- 
ences, see for examples Ref. jB| and Chapter 20 of Ref. Q . 

D. Superconducting Quantum Integrated Circuits 

There is no reason why the inductive loops cannot couple 
qubits that are far apart. In addition, a single qubit can 
be coupled to several other qubits as shown in Fig. |§|. This 
arrangement requires separate AC control lines for each of 
the qubits, which then demands localized on-chip oscillators. 
One can build up essentially arbitrary integrated circuits of 
superconducting qubits by this method. This flexibility in the 
construction of quantum computer architectures is one of the 
benefits of using superconducting Josephson junction circuits 
to perform quantum computation. The quantum integrated 
circuit could be set up to provide a number of useful features. 
For example Kj, one might be able to design the circuit and 
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A method for coupling a single qubit to other 



interactions in such a way that it automatically implements 
an intrinsically fault-tolerant quantum computer architecture 
such as those proposed by Kitaev ^] and Preskill ^lj. In 
addition, since the circuits are parallelizable in that different 
quantum logic operations can be performed in different places 
simultaneously, the circuit could be designed to provide the 
maximum possible parallelization of a particular problem such 
as factoring |Q , database searc h |r33]| , or computing a discrete 
quantum Fourier transform [p2|J 



VIII. DECOHERENCE 

We have shown how superconducting circuits can be used 
to construct qubits and quantum logic circuits. These super- 
conducting qubits have been idealized in the sense that we 
have ignored the effects of manufacturing variability, noise 
and decoherence. Manufacturing variability can be compen- 
sated for as discussed above: before performing any quantum 
computations, the properties of individual qubits can be mea- 
sured, recorded in a look-up table in a conventional computer, 
and used either to supply compensating calibration fields or 
to alter the frequencies with which control pulses are supplied 
to the qubits. 

From the point of view of the ultimate performance of a 
superconducting computer, a more pressing issue is that of 
environmentally induced noise and decoherence. In real sys- 
tems the performance of a qubit will be limited by dissipative 
mechanisms which cause the quantum state to decohere in 
time Td. The 'quality factor' for a qubit is the decoherence 
time divided by the amount of time it takes to perform fun- 
damental quantum logic operations The quality factor 
gives the number of quantum logic operations that can be 
performed before the computation decoheres, and should be 
10 4 or greater for the quantum computer to be able to per- 
form arbitrarily long quantum computations by the use of 
error-correction techniques p^-|69[|. 

Decoherence can be due to "internal" dissipation (quasi- 



particle resistance), or coupling to an environmental degree 
of freedom. It is also possible to couple to an environmental 
degree of freedom, without a dissipative mechanism, that will 
still lead to decoherence |70| 

We will now discuss some of the major sources of decoher- 
ence. 

Normal state quasi-particles can cause dissipation and en- 
ergy relaxation at finite temperatures in Josephson junctions. 
However, mesoscopic aluminum junctions have been shown 
to have the BCS temperature dependence for the density of 
quasi-particles. At low temperatures this density is exponen- 
tially small jnj, so quasi-particle tunneling will be strongly 
suppressed at low temperatures and at low voltages, as was 
seen in a system with multiple superconducting islands in 
Ref. We estimate a lower bound of 10 4 for the quality 

factor, given a sub-gap resistance of 10 10 Q, |7lj . 

The qubit can also decohere by emitting photons. We esti- 
mate this by classically estimating the rate photons are emit- 
ted by magnetic dipole radiation from oscillating current in 
the loop defining the qubit. For a loop with radius R with 
an alternating current of frequency v and rms amplitude of 
I m , the power transmitted to free space is P m — K m R 4 I'm^ 4 
where K m — 8/x 7r 5 /3c 3 . A typical rate for photon emis- 
sion is 1/tm = Pm/hv, which gives an estimate of the de- 
coherence time of t m = h/K m R i I r \v i . An estimate for the 
frequency is the Larmor frequency, (other characteristic fre- 
quencies such as the Rabi frequency are even smaller). For 
our qubit R ~ 1 /im, v w 10 GHz and I m = 1 nA, we find 
that t m ~ 10 7 s, so that this is not a serious source of deco- 
herence. However, it should be noted that proposals for using 
RF SQUIDs for qubits, involve currents of the order of 1 fiA 
and and loops of the order of 10 [im. These RF SQUIDs have 
t m ~ 10 -3 sec, which is substantially lower than for our qubits 
which can be made much smaller have much less current. 

Inhomogeneity in the magnetic flux distribution can also be 
a source of decoherence. This is similar to T2 in NMR systems. 
We estimate this for our system by calculating the amount of 
flux a 1 /im x 1 /im wire carrying 100 nA of current induces 
in a loop of the same size which has its center 3 /im away. 
We find that the induced frustration is about Sf = 10 -7 . 
If this is taken as an estimate of the typical variance of the 
frustration that difference qubits experience, then there will 
be a spread of operating frequencies among the loops. An 
estimate of td is the time for the extremes of this frequency 
differ by ir. This results in td ~ ir/(2ri8f), where we have 
taken the larger value from Eqn. B5l With r\/h ~ 600 GHz, 
we find td ~ 1.5 msec. The dipole-dipole interaction between 
qubits gives a time of the same order. 

We have also estimated the magnetic coupling between the 
dipole moment of the current loops and the magnetic mo- 
ments of the aluminum nuclei in the wire. At low tempera- 
tures where the quasi-particles are frozen out, the decoherence 
time for a single qubit is of the order of Ti which is exponen- 
tially large in the low-temperature superconducting state. For 
an ensemble of qubits, the decoherence time time may be of 
the order of milliseconds due to the different configurations of 
nuclear spins in the different qubits. However, this effect may 
be reduced by aligning the spins or by applying compensating 
pulse sequences. 

Coupling to ohmic dissipation in the environment has been 
modeled for superconducting qubits operating in the charging 
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regime |19J. They found that source of decoherence could be 
made sufficiently small such that the quality factor is large 
enough. Similar calculations for qubits in the superconduct- 
ing regime of circulating currents have not yet been done. 
Experiments to measure this decoherence time in our circuits 
are underway. In practice electromagnetic coupling to the 
normal state ground plane can limit coherence [ psf ; however, 
a superconducting ground plane can greatly reduce this cou- 
pling. 

Other possible sources of decoherence are the effects of the 
measuring circuit, the arrangement and stability of the con- 
trol lines for the magnetic fields, the ac dielectric losses in the 
substrate at microwave frequencies. These and other source 
of decoherence will have to be estimated in a real circuit en- 
vironment and measured. 

Taking 0.1 msec as a lower bound on the decoherence time 
and 10 nsec as a switching time, we find that the quality fac- 
tor is of the order of 10 4 . Furthermore, if the proper set of 
topological excitations is used to store information, the de- 
coherence time for quantum computation can be made sub- 
stantially longer than the minimum decoherence time for an 
individual junction circuit [BOj. 



IX. SUMMARY 

In this paper we have discussed a superconducting qubit 
which has circulating currents of opposite sign as its two logic 
states. The circuit consist of three nano-scale Josephson junc- 
tions connected in a superconducting loop and controlled by 
magnetic fields. One of the three junctions is a variable junc- 
tion made as a SQUID loop. This qubit has quantum states 
which are analogous to a particle with an anisotropic mass 
moving in an two-dimensional periodic potential. Numeri- 
cal calculations of the quantum states of the qubit have been 
made as well as physical estimates from a tight-binding ap- 
proximation. The advantages of this qubit is that it can be 
made insensitive to background charges in the substrate, the 
flux in the two states can be detected, and the states can be 
manipulated with magnetic fields. Coupled systems of qubits 
are also discussed as well as sources of decoherence. 
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APPENDIX A: CLASSICAL STABILITY 

In this appendix we find the eigen values of the stability 
matrix for the three junction potential and the range of frus- 
tration around / = 1/2 where there are two stable classical 
solutions with opposite circulating currents. 

The potential energy of the Josephson energy of the three 
junction qubit is given by Eqn. |l| 



U = = 2 + a — cos tpi — cos if2 — o. cos(27r/ + ip\ 
Ej 



¥"2) 
(40) 



We are interested in minimum energy phase configurations; 
that is, stable solutions of the following system of equations: 



dU 
dipi 
dU 
dip 2 



= sin tfii + a sin(27r/ + ip\ — if 2) = 

= sin if 2 — a sin(2-7r/ + ipi — if 2) — (41) 



The solutions {<f*,ifi2) complies with: sin<£* = — sin^>2 
sin if*. Then 



sin cp* — — a sin(2-7r/ + 2if*) 



(42) 



In order to check the character of the solution we compute 
the eigen values of the stability matrix, g ®.g v . , where 

d 2 u 



2 = cosi/Pi + acos(27r/ + <p\ — ^2) 
oif 1 

d 2 U ,„ , , 

— ? = cos (f2 + a cos(27r/ + ipi — if 2) 
dfi 



(43) 



d 2 U 
df\df>2 



-acos(27r/ + (pi — if 2) 



For the states with cos if* = cos if 2 = cos if* (these are the 
ones we are interested here), the eigen values are 



Ai = cos if* 

X2 — cos if* + 2a cos(2-7r/ + 2ip* ) 



(44) 



When / / 0, 1/2 we have used relaxation methods for com- 
puting if* . Both eigen values are greater than zero which as- 
sures the minimum energy condition. Fig. ^| shows the energy 
of the minimum energy configurations for a = 0.8. We find 
that there exists a region of values of the field for which two 
different minimum energy phase configuration coexist. 

Next we calculate the critical values of the external field for 
this coexistence. We can restrict our analysis to the region 
around / = 0.5; that is, [0.5 - ,f c ,0.5 + / c ] (where f e > 0). 
These extrema values of the field correspond to solutions for 
which one of the eigen values is positive and the other equals 
zero. The inset of Fig. |^ shows f c (oe). 

We first calculate f c when a > 1.0. The first eigen value 
which equals zero is Ai- Then at / = 0.5 ± f c , Ai = which 
implies if* = T n /2 mod 2n (here and below we associate the 
sign in f c with the sign of the phase in order to have f c > 0). 
Then, going to Eqn. WA we get 



sin(q=7r/2) = — asin(7r ± 2-nf c =p 71") 
±1 = ±asin(27r/ e ) 



(45) 
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and 



, 1 .1 
Jc = — arcsm — 
Itx a 



(46) 



We now calculate f c when 0.5 < a < 1.0. Now the first 
eigen value to equal zero is A2, and we have to solve: 

sinip* = — asin(27r/ + 2ip*) — asin(±27r/ c + 2ip*) 

cos(p* = -2acos(2yr/ + 2tp*) = 2a cos(±2tt/ c + 2ip") (47) 

We will use A = ±2irf + 2(p* , so that 

1 = sin 2 if* + cos 2 ip* 



= a 2 sin 2 A + 4a 2 cos 2 A 
= a 2 + 3a 2 cos 2 A 



(48) 



Then 

cos A 

COS if* = 



1 - a 2 



3a 2 



1-a 2 



A = =p arccos( 



= =p arccos(2 



1-a 2 



3q 2 



1 - a 2 



(49) 



Here we have followed the solution corresponding to cos(y?*) > 
0. Finally we have the solution for f c (A = ±2nf c + 2ip*) 



2tt 



2 arccos(2 



1-a 2 



arccos( 



1-a 2 . 
3a 2 



(50) 



APPENDIX B: TIGHT-BINDING ESTIMATE OF 
COEFFICIENTS OF CHANGE 

Recall that f a — f2 and fb~fi+ fil2. Assume that we 
change f a and fb independently. The minima in U occur at 
ifp — and tp* m = ±v?m Therefore, the energy due to the 
potential energy is for each of the minimum 



U_ 

Ej 



2 + 2/3 — 2 cos Lp* m — 2/3 cos(7r/ a ) cos(2ivf b + 2ip* m ) . 



(51) 



The change in the magnetic flux f a by Sf a causes a change in 
U of 



dU 

dfa 



5 f a = -27T/3 sin Tvf a cos 2ip° m Sf a 



(52) 



which is the same for the minimum at iip^i- Whereas, the 
flux fb causes a change 

dU 

TTT S fb = =F47T/3 COS TV f a Sill 2ip° m 8 f b 

dfb 

which has opposite signs for the two minimum. Therefore, 

— — = ~2n/3 sin nf a cos 2tp° m 5f a 1 — 47r/3 cos irf a sin 2f° m Sfb &z 
Ej 



(54) 

Recall that AF is the change is the energy between the two 
states when there is no tunneling. This is the second term in 



Eqn. |54| since the first term is only a constant for both levels, 
so that 



AF 

s7 



-4-7T/3 cos 7r/ a sin 2</5„ Sfb <j z 



(55) 



For this change AF = n5i + riSi; and since Sfb = 8\ + S2/2, 
we have n = 2r2 and 



— = 47T/3 cos ty fa sin 2ip° m 
Ej 



(56) 



We have found previously that costp° m = l/2a where a 
2/3 cos it f a so that with /„ = 1/3, 



= 2ny/l - 1/(4/3=) . 

Ej 



(57) 



To find the changes in At, we see that the changes in ti 
(huj m /2Tv)e- Sl/R are dominated by changes in Si, so that 



h ^ dfi 



(58) 



The changes in fb do not change Si to first order. Hence, 
changes in t come from changes in f a = fi only, so that 
81 — 0. But changes in f a are equivalent to changes in a in 
the three junction problem, so we can use Eqns. |5^ and |l6| 
and the fact that 2/3cos(7r/ a ) plays the role of a to find 



At = 5^(20 sin 7r/ a )*/„. 



(59) 



This allows us to write S2 = rft^/Ej /E c where r\ is of the 
order of unity. For the operating point we find r\ ~ 3.5. There- 
fore, changes in H due to changes in ti go like cr x . These 
tight-binding estimates for /3 = 0.8 and f a = 1/3 give si = 
and S2 — 0.03. 
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